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Abstract. In this paper we use the pseudofermion dynamical theory (PDT) in 
the study of the one-electron removal singular spectral features the one-dimensional 
Hubbard model. The PDT reveals that in the whole (k, w)-plane such features are 
of power-law type and correspond to well defined lines of three types: the charge 
singular branch lines, the spin singular branch lines, and the border lines. One of 
our goals is the study of the momentum and energy dependence of the spectral- 
weight distribution in the vicinity of such lines. We find that the charge and 
spin branch lines correspond to the main tetracyanoquinodimethane (TCNQ) peak 
dispersions observed by angle-resolved photoelectron spectroscopy in the quasi- ID 
organic conductor tetrathiafulvalene-tetracyanoquinodimethane (TTF-TCNQ). Our 
expressions refer to all values of the electronic density and on-site repulsion U. The 
weight distribution in the vicinity of the singular spectral lines is fully controlled by the 
overall pseudofermion phase shifts. Moreover, the shape of these lines is determined 
by the bare-momentum dependence of the pseudofermion energy dispersions. 



PACS numbers: 71.20.-b, 71.10.Pm, 72.15.Nj, 71.27.+a 
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1. INTRODUCTION 

Early studies of quasi-one-dimensional (ID) compounds have focused on the various 
low-energy phases, which are not metallic and correspond to broken-symmetry states 
PU El- Recently, the resolution of photoemission experiments has improved, and 
the normal state of these compounds was found to display exotic spectral properties 
jH El • The study of the microscopic mechanisms behind such properties remains until 
now an interesting open problem. Indeed, the finite-energy spectral dispersions recently 
observed in such metals by angle-resolved photoelectron spectroscopy (ARPES) reveal 
significant discrepancies from the conventional band-structure description [?,C1II1E1III]- 
There is some evidence that the ID Hubbard model pj E] describes successfully 
the transport properties and other exotic properties observed in some low-dimensional 
materials [0] and that the electronic correlation effects described by the model could 
contain the finite-energy microscopic mechanisms [3J E] that control the above finite- 
energy spectral properties. Until recently very little was known about the finite-energy 
spectral properties of that model for finite values of the on-site repulsion U. This is in 
contrast to simpler models [TU]. Indeed, usual techniques such as bosonization jTTj and 
conformal-field theory [T21 EEH1 CI EES EEE1 EE3 CHI EH H3 do not apply at finite energy. 
Valuable qualitative information can be obtained for U — > oo by use of the method 
of Refs. [22l EHj. However, a quantitative description of the finite-energy spectral 
properties of quasi-lD metals requires the solution of the problem for finite values of 
the on-site Coulombian repulsion U. The method of Ref. [21] refers to features of the 
insulator phase. For U ~ 4t, where t is the transfer integral, there are numerical results 
for the one-electron spectral function [23j. Unfortunately, the latter results provide very 
little information about the microscopic mechanisms behind the finite-energy spectral 
properties. 

Recently, the preliminary use of the finite-energy holon and spinon representation 
introduced in Ref. PS] and of the related pseudofermion description of Refs. 
|27[ 12*5] . revealed that most singular features of the one-electron removal spectral 
function correspond to separate charge and spin branch lines jH Ej- Interestingly, 
the singular branch lines associated with the one-electron removal spectral function 
show quantitative agreement with the TCNQ peak dispersions observed by ARPES in 
the quasi-lD organic conductor TTF-TCNQ [HE1E]- However, these studies provide 
no information about the momentum and energy dependence of the spectral-weight 
distribution in the vicinity of the charge and spin branch lines. A preliminary study of 
that dependence was recently presented in short form in Ref. [23 • Shortly after a study 
of the same problem by means of the dynamical density matrix renormalization group 
(DDMRG) method led to very similar results [3U] . 

The main goal of this paper is the extension of the preliminary results presented 
in Ref. [211] • Our investigation relies on the pseudofermion dynamical theory (PDT) 
introduced in Ref. j2Z]- Such a theory provides the general finite-energy spectral- 
weight distributions for the metallic phase of the model for all values of energy and 
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momentum. However, the amount of one-electron removal weight in regions away from 
the singular spectral lines is small. Therefore, here we limit our studies to the vicinity of 
such lines. Our investigation focus on the one-electron removal spectral function for all 
values of U and electronic density in the vicinity of the singular charge and spin branch 
lines, which turn out to be the most important spectral features for the description of 
the unusual TCNQ photoemission spectral lines. The evaluation of the small spectral- 
weight distributions away from the singular spectral features considered here requires 
the use of involved numerical calculations which will be carried out elsewhere. 

The pseudofermion description refers to the pseudofermion subspace (PS) j^Zj- The 
one-electron removal excitations studied in this paper are contained in such a Hilbert 
subspace. The PDT refers to all energy scales but in the limit of low energy leads to the 
known conformal- field-theory spectral-function and correlated-function expressions |31j . 
In Ref. the PDT is combined with the renormalization group for the study of the 
microscopic mechanisms behind the phase diagram observed in the (TMTTF) 2 X and 
(TMTSF) 2 X series of quasi-lD organic compounds. These studies are consistent with 
the latter phase diagram and explain why there are no superconducting phases in TTF- 
TCNQ. The PDT used in this paper in the study of the one-electron removal spectral 
function of the ID Hubbard model applies to related integrable interacting problems 
[3H] and therefore has wide applicability. 

The paper is organized as follows: In Sec. I we introduce the one-electron removal 
spectral-function problem, the ID Hubbard model, and some basic information about 
the pseudofermion description. The introduction to the PDT and of the corresponding 
general spectral-function expressions used in our investigation is the subject of Sec. 
II. Section III is devoted to the the study of the weight distributions for one-electron 
removal. Moreover, we investigate the limiting behavior of the spectral function in the 
vicinity of the branch lines both for U/t — > and U/t — > oo. The discussion of the 
relation of our theoretical predictions to the TCNQ branch lines observed by angle- 
resolved photoelectron spectroscopy in the quasi- ID organic compound TTF-TCNQ 
and the concluding remarks are presented in Sec. IV. 

1.1. THE MODEL AND THE ONE-ELECTRON REMOVAL SPECTRAL 
FUNCTIONS 

The ID Hubbard model reads, 

H = -tJ2[ci*Cj+i,<r + h.c] + U^2h^h jti , (1) 

where c] ^ and c 3j(T are spin-projection a =f, j electron operators at site j = 1, 2, N a 
and hj )(J = q a Cj ia . The model (JTJ) describes spin- up electrons and N± spin-down 
electrons in a chain of N a sites. We denote the electronic number by N = + N[. 
The number of lattice sites N a is even and very large. For simplicity, we use units such 
that both the lattice spacing and the Planck constant are one. In these units the chain 
length L is such that L = N a . Our results refer to periodic boundary conditions. We 
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consider an electronic density n = n-j- + nj_ in the range < n < 1 and a spin density 
m = n-f— nj_ such that m — > 0, where n CT = N a /L and a =|, I- However, the calculations 
are performed for finite values of the spin density m. For U/t > the limit m — > is 
taken in the end of the calculation and leads to the correct m = results. The Fermi 
momentum is kp = irn/2 and the electronic charge reads — e. 

The one-electron removal spectral function B(k, uj) is given by, 

B(k,uj)= \(f\c kl<r \GS)\ 2 5(u + E f -E GS ), lu<0. (2) 

CT / 

Here Cfc jCr is an electron annihilation operator of momentum k and \GS) denotes the 
initial iV-electron ground state. The / summation runs over the N — 1-electron excited 
energy eigenstates and [Ef — Eqs] are the corresponding excitation energies. We use an 
extended momentum scheme such that k G (— oo, +oo) for the expression given in Eq. 
(J2J), yet it is a simple exercise to obtain the corresponding spectral function expression 
for the first Brillouin zone. 

1.2. THE PSEUDOFERMION DESCRIPTION 

The pseudofermion description [2E1 E3 EHl EH] involves cv pseudofermion branches such 
that v = 0,1,2,... and sv pseudofermion branches where v = 1,2,.... However, most 
of our study refers to the two branches that contribute to the one-electron removal 
dominant processes considered below. For the one-electron removal spectral function 
corresponding to densities < n < 1 and m — > the above-mentioned dominant 
processes involve the cO pseudofermions and si pseudofermions of Ref. [2Zj. The 
pseudofermions are related to the pseudoparticles studied in Refs. E3] by a 

canonical transformation. Such a transformation introduces shifts of order 1/L in the 
pseudoparticle discrete momentum values and leaves all other pseudoparticle properties 
invariant. As a result of such momentum shifts and in contrast to the pseudoparticles, 
the corresponding pseudofermions have no residual-interaction energy terms. 

Most of our final expressions refer to densities < n < 1 and m —>■ 0. According to 
the studies of Refs. |2E1 EH EH EE] the cO pseudofermion is a spin- less and 77-spin-less 
object that carries charge — e and the si pseudofermion is a charge-less and spin-zero 
two-spinon composite object. (That charge value corresponds to the description of 
the transport of charge in terms of electrons j2Ej-) An important point that follows 
from the above analysis is that while the cO pseudofermion branch is associated with 
excitations of charge character, the si pseudofermion branch describes spin excitations. 
The oiv = cO, si pseudofermions carry canonical-momentum q = q + Q® u (q)/L. Here q 
is the bare-momentum and Q^ u (q)/2 is the overall scattering phase shift given by |28j . 

Qt(q)/2 = n Yl J2®«u, a iu>(q,q')AN a , v ,(q'); av = c0,sl.(3) 

a'v'=cO, si q' 

On the right-hand side of this equation AN au (q) = AJ\f a v(q) is the bare-momentum 
distribution function deviation AN au (q) = N au (q) — N® u (q) and the elementary two- 
pseudofermion phase shifts ®av,a r v'(q, Q 1 ) in units of 71 are such that +7T^ aVja r l/ >(q,q') 
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(and —Tr$au,a'v'(q,q')) gives the phase shift acquired by the bare-momentum q av 
pseudofermion or hole wave function when such an object is scattered by a bare- 
momentum q' aV pseudofermion (and oiv' pseudofermion hole) created under a ground- 
state - excited-energy-eigenstate transition [28 . The bare-momentum distribution- 
function deviations AN a / u > (q') of Eq. (J3J) result from such a transition. In expression (J3J) 
these deviations refer to the cO and si branches only and thus describe excited energy 
eigenstates generated by one-electron removal dominant processes. The corresponding 
general expression is given in Eq. (2) of Ref. |2Hj and involves summations over all 
pseudofermion branches. 

Note that for each av branch the continuum bare-momentum q corresponds to a set 
of discrete bare-momentum values qj such that qj + i — qj = 2n/L. Here j = 1, 2, N* u 
and the number N* v = N au + N% v is given in Eqs. (B6)-(B8) and (Bll) of Ref. [H|. 
N av and N^ u denote the number of av pseudofermions and av pseudofermion holes, 
respectively. N* u equals the number of sites of the effective av lattice [2E], which 
plays an important role in the pseudoparticle and pseudofermion descriptions. For the 
av = cO, si branches the number N^ v reads, 



where in our case N^ x = N c0 — 2N S \. We used Eq. (B.ll) of Ref. to derive the 
expression given in Eq. (J3J for N^, such that Nf^ = for the m —>■ initial ground 
state. 

Although the av pseudoparticles carry bare-momentum q [2H1 EH], one can also 
label the corresponding av pseudofermions by that bare-momentum. Indeed, the latter 
pseudofermions carry canonical-momentum q = q + Q® u (q)/L, but this latter expression 
provides an one-to-one relation between the bare-momentum q and the canonical- 
momentum q. The pseudoparticles have residual-interaction energy terms which do 
not allow the expression of the electronic spectral functions as a convolution of av 
pseudoparticle spectral functions 27 . A property which plays a central role in the 
PDT is that for the corresponding av pseudofermions, such a residual-interaction energy 
terms are exactly canceled by the overall scattering phase shift Q* l/ (g)/2. By canceling 
the residual interactions exactly, the associated canonical-momentum shift Q^ v (q)/L 
transfers the information recorded in these interactions over to the pseudofermion 
canonical-momentum. It is found in Ref. that the av pseudofermion or av 

pseudofermion hole overall phase shift acquired under a ground-state - excited-energy- 
eigenstate transition is given by Q au (q)/2 = Q° au /2 + Q® u (q)/2. Here Q® u (q)/2 is 
provided in Eq. © and Q° av /2 = 0, ±tt/2 is a av pseudoparticle or hole overall 
scattering-less phase shift whose value is well defined for each excitation subspace 
spanned by energy eigenstates with the same pseudofermion numbers [23 EH] • It is such 
that under a ground-state - excited-energy-eigenstate transition the av pseudoparticle 
and hole discrete bare- momentum value qj is shifted by Q° av (qj)/L. The overall phase 
shift Qau(q)/2 leads to a corresponding canonical- momentum shift Q a u{qj)/L for the 
discrete canonical-momentum values of the av pseudofermions and holes. According 



oo 



= [N a - N c0 ] ; N* =N c0 -2Y, , 



(4) 
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to the studies of Ref. j2H], the overall phase shift Q av {q)/2 is associated with the olv 
pseudofermion or hole S matrix exp {iQ au (q)}. 

The pseudofermions have energy bands e c0 (g) and e s i(g) such that | q\ < n and 
| q\ < hp, respectively. These energy dispersions are plotted in Figs. 6 and 7 of Ref. 
[33] and are defined by Eqs. (C.15) and (C.16) of Ref. [2H], respectively. Also the group 
velocity v au (q) = de ai/ (q)/dq and the Fermi-point velocity v au = v au (q Fau ) appear in the 
spectral-function expressions. Here q Fc0 = 2k F and q Fsl — k F ^ define the ground-state 
Fermi points |2fi[ I27j. In Sec. Ill we confirm that the pseudofermion energy bands fully 
determine the shape of the one-electron removal spectral-function in the proximity of 
the branch lines studied in that section. In the ground state the si pseudofermion band 
is filled and the cO pseudofermions occupy the bare-momentum domain < | q\ < 2k F 
(leaving 2k F < | q\ < it empty). 

2. THE GENERAL SPECTRAL-FUNCTION EXPRESSIONS USED IN 
OUR STUDY 

The one-electron removal problem studied in this paper involves only the PDT dominant 
processes defined in Ref. [2Z]- The dominant (and non-dominant) processes correspond 
to the i = (and i > 0) terms on the right-hand side of the general spectral-function 
expression given in Eq. (41) of the first paper of that reference. Furthermore, only 
such dominant processes contribute to the one-electron removal spectral-function power- 
law expressions obtained in this paper. The initial ground state and excited energy 
eigenstates of Eq. (J2J) that are associated with the one-electron removal dominant 
processes can be expressed in terms of occupancy configurations of the cO and si 
pseudofermions. The first step of the evaluation of the spectral function (J2J) involves the 
expression of the electronic operators in terms of these quantum objects, as described 
in Ref. [27J. For one-electron removal, the dominant processes involve creation of one 
pseudofermion hole both in the bands e c o(g) and e s i(g). The excited energy eigenstates 
generated by these one-electron removal dominant processes belong to subspaces whose 
cO pseudofermion and si pseudofermion hole number ground-state deviations are given 
by, 

AiV* = -AiV c0 = 1 ; AA£ = 1 . (5) 

Moreover, as discussed in Sees. Ill and IV, the main one-electron removal spectral- 
function singular features are associated with the charge cO and spin si pseudofermion 
branch lines. 

The domain of the (k, c<j)-plane whose spectral weight is generated by one-electron 
removal dominant processes is contained in the region of negative uo/t values of Fig. 1 
of Ref. [22] • (That figure uses the extended momentum scheme also used here.) Such a 
domain is limited above by the s line for momentum values between k = and k = k F , 
c" line from k = k F until that line reaches the c' line at k = 2k F , d line from the 
momentum k = 2k F until k = 3&f, and s line between the momentum values k = 3kF 
and k = 5k F . The same domain is limited below by the lowest line of the figure. 
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Use of the general PDT distributions studied in Ref. for the one-electron 

removal problem involves all dominant processes associated with the pseudofermion- 
number deviations given in Eq. (|5]l. In addition to small- momentum and low-energy 
cO and si pseudofermion particle-hole processes, which conserve the pseudofermion 
numbers, such dominant processes involve creation of one cO pseudofermion hole and one 
si pseudofermion hole at bare-momentum values q and q', respectively. When both the 
bare-momentum values q and q' of the two created pseudofermion holes are away from 
the Fermi points and such that v c o(q) 7^ v s i(q'), the corresponding dominant processes 
do not lead to singular spectral features and generate the one-electron removal spectral 
weight for (k, u;)-plane regions away from these features. Since the latter weight is 
small, in this paper we do not study its intensity distribution in the (k, u;)-plane. A 
first type of singular feature corresponds to lines generated by such processes where both 
created objects move with the same group velocity, v c o(q) = v s i(q'), and the spectral 
feature corresponds to a border line, uj = UBhik) = [±e c0 (g) — ei s (q')] ^v c0 (q),v 3 i(g'), m the 
(k, o>)-plane. In this case the spectral function reads, 

B{k, u) « C BL {k) {u - uj BL (k)Y 1/2 , (6) 

in the vicinity and just above such a line. However, as discussed in Sec. IV, the only 
existing border line for the density of TTF-TCNQ leads to a weak TCNQ spectral 
feature mentioned in Sec. IV, which we do not study in Sec. III. 

The second type of spectral feature corresponds to the branch lines studied in this 
paper, such that either q or q' equals one of the corresponding pseudofermion branch 
Fermi points. Thus, such features are generated by processes where a <xu pseudofermion 
hole is created for all the available values of bare-momentum q and a second a'u' 
pseudofermion hole is created at one of its two Fermi points ± q% a i u i, where au = cO, si 
and a'u' = si, cO, respectively. According to the general theory of Ref. [21], f° r densities 
< n < 1 and < m < n the functional, 

^ = ajv2£ + 9A^hAy = (, A < t + 9iMA y , (7) 

where 1 = ±1 and au = cO, si controls the exponents of the spectral-function power- 
law expressions in the vicinity of the charge and spin branch lines. Here AiV^, L is 
the deviation in the number of au = cO, si pseudofermions at the Fermi points and 

Q*( L QFav)/2 an d Qav(t<qFav)/2 = Qav/^ + ( L QFav) / ^ are tne overall scattering 
and overall, respectively, phase shifts of au pseudofermion or hole scatterers at the 
bare-momentum Fermi values iq% av = ±qp aiJ . Note that the deviation AN^ U L = 
lAN^ p l + Q° au /2ir involves both contributions from the scattering- less overall phase 
shift and the number deviation AN®' I f l generated by the creation or annihilation 

pseudofermion processes at the Fermi points. We also consider the au = cO, si current 
number deviation 2AJ% U = AN£ V>+1 - AN^_ V 

For excited energy eigenstates generated by processes involving pseudofermion 
occupancies in the vicinity of the cO or si Fermi points, the functional (JJJ) and the 
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corresponding branch-line exponent expressions involve the following Fermi-point two- 
pseudofermion phase-shift parameters, 

£LaV = <W ^,v> + H W) ®<xv,a'v>(qFav,t<<lFa>v>) I J =0, 1 , (8) 

i=±l 

where ctz/ = cO, si and aV = cO, si. In the limit m — > 0, these parameters are given by 

£cOeO = VCoj £cOsl = 0' £slcO = "V^i £slsl = V^; £cOcO = Co> £cOsl = Co/2, CslcO = 0' 

and Csisi = 1/V2- Here £o is the parameter defined in Eq. (74) of Ref. [?] and in 
the text above that equation. It is such that £o ~~ > \/2 and £o — > 1 as Z7/i — > and 
C//t — > oo, respectively. 

We emphasize that the limits m — > 0, U/t —*■ and Z7/t — > 0, m — > Odo 
not commute and lead to different values for the parameters (JBJ): while for m — ► 0, 
17/* - one finds £ c ° 0c0 - 1/^2, & sl - 0, £° lc0 - "Vv 7 ^, £ s ° lsl - A Cko - v 7 ^, 
^cosi V V2, ^lico 0, and £ lal -> l/v 7 ^, for C//t -> 0, m -> the result is £? 0c0 -> 1, 
£cOsi — *■ 0) £sico — *■ 0' an d £sisi ~^ 1 both for j = and j = 1. Also the phase shifts 
®si,co(kFi, q) and <& s i,si(&f.|.j 9) have different values in the limits m — > 0, £7/7 — > 
and — > 0, m — > 0, respectively. Therefore, for some excitations the spin functionals 
2Aji and 2Aj L 1 have different values in these two limits. For such excitations we provide 
below the spectral-function exponent associated with the limit U/t — *■ 0, m — > 0, which 
is that which leads to the correct U/t = spectral-function behavior for m — 0. This 
can be confirmed by studying the limit of the corresponding exponents as U/t — > for 
m > 0. 



£1. GENERAL EXPRESSIONS IN THE VICINITY OF THE PSEUDOFERMION 
BRANCH LINES 

Here and in section III we label the olv pseudofermions by their bare-momentum q, 
which is related to the canonical-momentum as q = q + Q^ v (q)/L. The momentum 
values k of the one-electron removal olv branch lines points (k, luj av {k)) are determined 
through the bare-momentum value q of the created civ pseudofermion or hole by the 
following parametric equations, 



k 
h 



k -q 



Q 



k + k 



-600/(9) 



4k F AJ* + 2k Fl AJ? 1 . 



(9) 



According to the studies of Ref. [2Zj , the one-electron removal spectral function (J2J) 
has in the vicinity of the olv branch lines the following power-law expression for finite 
values of the energy uj < such that —(u+uj au (q)) is small and positive and u av (q) 7^ 
where the point (k, —uj av {q)) belongs to the branch line, 

[uj + UJ au (q)]\( av (q) 



B(k } uj) « C av {q) 



= 8{u> + u av {q)}] ( au (q) = -l; av = cO, si . (10) 

Here the first expression is the leading-order term of a power-law expansion in the small 
energy deviation — [uj+uj av {k)\ relative to the branch-line energy whose exponent ( au (k) 



The Hubbard model description of the. 



9 



and pre- factor C au (k) have the following general form, 

Cau(q) = - l + Co(?) ; 

Co(q) = 2A+ 1 (g) + 2A- 1 (g) + 2A+ 1 (g) + 2A; 1 1 (g) > ; 




Vav(q)\ 
Vau{q)\ 

i a u{q) 



■ (11) 



In these expressions I au (q) = [l+[u)+u) au (q)}/Q}, the small positive energy Q corresponds 
to the energy range of the small-momentum and low-energy cO and si pseudofermion 
particle- hole elementary processes, Fq(z) is an even function of z which is defined in 
Ref. [27|, and the functional 2A L av is given in Eq. ( |7| ). (In the numerical evaluation 
of the spectral-function expressions the value of Q is determined by the normalization 
procedure associated with imposing the k and u spectral-function sum rules j^Zj; The 
expression provided in Ref. |2Z1 for F (z) is valid when the four parameters 2A^ are 
finite; An expression for F (z) valid when some of these parameters vanish will be given 
elsewhere.) A full quantitative study of the pre-factor C au (q) whose general expression 
is provided in Eq. involves the numerical derivation of the lowest-peak weight 

functional introduced in Ref. j2ZJ, which is included in the expression of the function 
Fq(z). Such a quantitative study, which requires involved numerical calculations, is 
beyond the goals of this paper and will be carried out elsewhere. However, the general 
C a u(q) expression of Eq. (fTTj) is useful for our studies, once it can be used to extract 
information about the behavior of the pre-factor C au (q) as U/t — > and find out whether 
for finite values of U/ 1 that function is vanishing or finite and also what its relative value 
for different branch lines is. 

A av branch line whose exponent Cav{q) is negative for a given domain of k = 
— [ko — q] values is called a singular branch line. In this case the weight distribution 
shows a singular behavior at the branch line, and we expect that the spectral peaks will 
be observed in a real experiment. This was confirmed for the present case of one-electron 
removal, as discussed in Sec. V. On the other hand, when for a (k, w)-plane region in 
the vicinity of the branch line and contained inside the one-electron removal dominant- 
weight domain of Fig. 1 of Ref. [29J the exponent (fTTj) is such that < ( a v(q) < 1) 
the spectral feature refers to an edge branch line. Finally, < Cav(q) < 1 for regions 
away from that domain and ( a v{q) > 1 for any (k, u;)-plane region are in general a 
sign of near absence of spectral weight. For one-electron removal and thus uj < 0, the 
singular and edge branch lines are represented in Fig. 1 of Ref. |2E] by solid and dashed 
lines, respectively. The dashed-dotted lines of that figure are either limiting lines for the 
domain of weight generated by dominant processes or lines associated with exponents 
larger than one. The lowest limiting line of the figure corresponds to a singular but 
weak spectral feature called a border line, as mentioned in Sec. V. 

The general branch-line spectral-function expressions defined by Eqs. (jlOjl and (jllj) 
are not valid in the vicinity of the low-energy branch-line end points. These end points 
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correspond to values of k and u such that to ~ iv av {k + k$ ) where, 
ko = k - l q° Fav ; cm/ = cO, si , 6 = ±1 , 



and fco is the momentum given in Eq. ©• In this low-energy limit the physics is that of 
the so called low-energy Tomanaga-Luttinger liquid (TLL) regime, where bosonization 
[HI and conformal-field theory [1211311111131131111111111111201 are applicable. 



2.2. LIMITING LOW-ENERGY BEHAVIOR NEAR THE BRANCH-LINE END 



The av = cO, si branch lines also exist for small positive values of — u. Such a regimen 
corresponds to the above-mentioned values of k and u such that uj « iv av {k + k$). In 
this case the expression (|10p does not apply and instead the general PDT leads to the 
following one-electron removal spectral- function expression |27j . 



where av = cO, si, l = ±1, and the pre-factor C av>L is a real and positive number. 
The low-energy spectral function expression given in Eq. (|12|) refers in general to the 
proximity of a av — cO, si branch-line end point. In turn, the finite-energy expression 
(IK)j) applies in the vicinity of the av = cO, si branch lines when the av = cO, si branch- 
line group velocity v au (q) is such that v au (q) ^ iv av . As the value of v au (q) approaches 
that of the Fermi point velocities ±v au , v ai/ (q) — > LV au , the spectral function k and uj 
values reach the vicinity of a branch-line end point and thus the spectral function is 
given by expression (|T2j) instead of (JTUJ) • 

Interestingly, for uj w iv av [k — lk^) in the momentum expression ©, the validity 
of the spectral- function expression given in Eq. ()12j) corresponds to the TLL regime. 
Therefore, in this limit the above exponent — 1 + C^ u L provided by the general PDT of 
Ref. [2Z] must equal that corresponding to the low-energy universal TLL expressions. 
Indeed, it is straightforward to show that the general exponent — 1 + (® u L of expression 
(JT2J) is identical to that given in Eq. (5.7) of Ref. [I2j. When applied to specific spectral 
functions such that u < 0, expression (JT2J) provides the universal and well known low- 
energy TLL behavior for the ID Hubbard model [I21[I3[Tl[I3[T3[TIl[T3[T3[lIj, 
Tomonaga-Luttinger model [?, ?, ?], and many other models whose low-energy physics 
corresponds to the same universality class. When — 1 + £° < 0, such an expression 
refers to a linear singular spectral feature. 

There is a cross-over region between the finite-energy and low-energy regimens 
corresponding to the spectral-function expressions (fTUj) and (fT2j) . respectively. The 
momentum and energy width corresponding to such a crossover regimen is very small 
and is fully controlled by the value of \v av (q) — iv av \. The low-energy TLL behavior 
emerges when \v av (q) - iv au \ « \a a u(q°F a v)(k + fcjf)|, where a av (q) = dv av (q)/dq for 
av = cO, si. As the value of the branch-line bare-momentum q of Eq. (JUJ) approaches 
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l Qfoo> the behavior (|T2j) is reached. Importantly, for smaller values of \a ai> (q Fau )\ the 
value of \v av (q) — tv av \ can remain small for larger values of \k + k$\ and thus of 
uj ~ iv av ik + fcjf). It follows that the momentum and energy widths of the (k,u)- 
plane region in the vicinity of the point (— , 0) where the TLL liquid behavior is valid 
increase for decreasing values of \a a u(Qp au )\, provided that v av is finite. For instance, in 
the limit of zero spin density, m — > 0, the value of |a s i(g)| is small in two relatively large 
q regions in the vicinity of q = —k F and q = +k F , respectively, and thus the domain of 
the corresponding spin si branch lines where the TLL expression ()12j) is valid increases 
in that limit. (We recall that the general exponent (JTT|) cannot be obtained by the TLL 
low-energy methods.) 

From the use of conformal- field theory [?], it is known that when the above branch- 
line end point (— , 0) is reached through low-energy lines other than the above branch 
lines the spectral-function expression is different from 1)1 2j) . Examples of such points 
are the points (k = kp, u = 0) and (k = 3k F , uj = 0) shown in Fig. 1 of Ref. The 
general theory provides spectral-function expressions that are valid when these points 
are approached by lines that are contained in the finite-weight regions and do not cross 
the av branch lines ending at the same points. The velocity v = 1/z = (u/ (k + k^) plays 
an important role in these expressions. Indeed, following the studies of Ref. |2Z|, in the 
proximity of the point (— k^ , 0) but for values of k and uj such that uj fa v (k + k^), where 
v 7^ ±f s i, ±^co and depending on the specific point, the z — 1/v domain is bounded 
by two of the four values —l/v c0 , —l/v s i, l/f s i> an d l/^cQ, the momentum and energy 
weight-distribution dependence has the following general expression for finite values of 



where F (z) is the function appearing in Eq. ((TTJ) , which is defined in Ref. j2Z], and (q 
is the above functional Co = 2A+Q 1 + 2A ( T 1 + 2A+! 1 + 2A;/. Here 2A^ is the functional 
of Eq. (JZJ) for the bare-momentum distribution function deviations associated with 
the excitations that control the spectral-weight distribution in the vicinity of the point 
(— fcjf, 0). We note that the u dependence of the spectral-function ()T3)l is that also 
obtained by conformal- field theory in Ref. [?] and by another method in Ref. [?]. 

Since the power-law spectral-function expressions (fi^j) and (fT^j) were studied 
previously by other methods and refer to the vicinity of low-energy lines and isolated 
zero-energy points in the (k, c<j)-plane, respectively, that are not of interest for the low- 
energy phase of the organic compound TTF-TCNQ, here we limit our study to the 
more complex problem of the finite-energy branch- line spectral weight. Indeed, the 
low-energy phase of TTF-TCNQ is not metallic and refers to broken-symmetry states 
[H El Ej whereas the branch-line spectral features given by Eq. (jl())l refer to finite-energy 
values which correspond to the unusual metallic state of that organic compound. 
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3. THE ONE-ELECTRON REMOVAL BRANCH LINES 

In this section we use the general branch-line expressions and associated quantities 
considered above in the study of the one-electron removal finite-energy spectral-function 
singular and edge branch lines. The ground-state - excited-energy-eigenstate transitions 
to the subspace whose pseudofermion number deviations are given in Eq. (j5J), generate 
several au pseudofermion branch lines whose location in the (k, u;)-plane is shown in 
the Fig. 1 of Ref. for uj < 0. The s = si branch-line shown in the figure, which 
connects the points (k = —kp, uj = 0) and (k = fcp, uj = 0), is generated by creating 
the cO pseudofermion hole at one of its Fermi points, and the si pseudofermion hole 
for bare-momentum values in the domain defined by the inequality \q\ < kp. We 
emphasize that in addition to creation of a cO pseudofermion hole at q = 2k p (and 
q = —2kp), this excitation includes a collective bare-momentum shift Q%/L = +7i/L 
(and Qco/L = —n/L) for the whole cO pseudofermion Fermi sea. 

By considering the same processes, plus transferring a cO pseudofermion from the 
Fermi point —i'2kp to the Fermi point t'2kp, two other si branch lines are generated, 
which connect the points (k = —i'3kp, uj = 0) and (k = —L f 5kp, uj = 0) where l' = ±1. 
The l' = —1 line is labeled by s in Fig. 1 of Ref. [2H], where it appears for uj < 0. On 
the other hand, there are four cO pseudofermion branches lines which connect the points 
[k = —3kp, uj = 0) and {k = kp, uj = 0), {k = —kp, uj = 0) and {k = 3kp, uj = 0), 
{k = —5kp, uj = 0) and {k = —kp, uj = 0), and (k = kp, uj = 0) and {k = 5kp, uj = 0). 
The first, second, and fourth of these lines are labeled by c, c', and c", respectively, in 
Fig. 1 of Ref. [29], where their k > parts are shown for uj < 0. Below we study the 
spectral-weight distribution in the vicinity of these seven one-electron removal branch 
lines. 

We start by evaluating the weight distribution corresponding to the first si 
pseudofermion branch line mentioned above. The specific form of the general expressions 
(JHJ) for the points (k, —uj s i(k)) belonging to the si pseudofermion branch line in the 
m —>■ limit, corresponds to ko = and reads, 

q = k; -uj s i(q) = e sl (q) . (14) 

Here e s i(q) is the energy dispersion given in Eq. (C.16) of Ref. |2E] and plotted for 
m — > in Fig. 7 of Ref. |34j . (The energy dispersions e c o(q) and e s i(q) appearing in 
other expressions of this section are those defined in Eqs. (C.15) and (C.16) of Ref. 
[2H] and plotted in Figs. 6 and 7 of Ref. [31], respectively.) We recall that the k > 
part of this si pseudofermion singular branch line is labeled by s in Fig. 1 of Ref. [29], 
where it connects the points (k — 0, u — e s i(0)) and {k = k F , uj = 0). In this case the 
general spectral-function expression (fTU|) applies provided that the specific expression 
associated with the excitations around the point (k, —uj s i(k)) of the functional 2A L au 
defined in Eq. ((7j) is used. This expression is a function of k = q and corresponds to 
the m —>■ limit of the following quantity, 

2A^( g ) = {-i, ^ - <WNL„ q)} 2 ; <*v = c o, si . (15) 
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Here the value of the two-pseudofermion phase shift <& a u,si an d that of the two- 
pseudofermion phase shifts & a u,a'v' appearing in other expressions of this section is 
uniquely defined in terms of the solution of a system of integral equations [2Z1 and the 
parameter £°„ c0 general expression is provided in Eq. (jHJ. 

Direct use of expression (|1U|) in the m — > limit, leads to the following expression 
for the one-electron removal spectral function, 

B(k,uj)^C sl {q){— j ; ( sl (q > -1 

= 6(u + uj sl (q)); C»i(?)=-1, (16) 

which corresponds to energy values just below the branch line for ( s i{q) > — 1 and at 
that line for ( s i{q) = — 1 arid to bare- momentum and momentum values in the range 
— kp < q < kp and — kp < k < kp, respectively The pre-factor C s i(q) given in Eq. 
(111)1 is finite for all values of the q domain, except in the vicinity of the branch line 
end points, where v s i(q) ~ ±i> s i and provided that Csi{o) > — 1 the spectral function is 
instead of the form given in Eq. (|12j). When ( s i(q) = — 1 the second spectral-function 
expression of Eq. (J16j) applies. It refers to the whole branch-line momentum domain. 

As the spin density m approaches zero, we find the following exponent expression 
valid for all values of U/t and electronic density n, 

U?)=-i+ E E + 

au=cO, sl t=±l 

= " 1+ E {^r + L^sii^kp, q )} 2 
L =±i 

+ J2{-7r7K + ^si,si(ikp,q)} 2 , (17) 

where the second expression was obtained by taking the limit m — ► in the first- 
expression quantities and the parameter £ is defined in Eq. (74) of Ref. [?]. The 
dependence of the exponent (j!7)l on the momentum k is obtained by combining Eqs. 
(114)) and (|17|) . The exponent £si of Eq. (|17|) is negative for all values of momentum 
and is plotted in Fig. 1 as a function of the momentum k for k > 0, several values of 
U/t, and electronic density n = 0.59. So the corresponding spectral-function expression 
(116)) describes a singular branch line. The exponent ( s of the figure is the exponent (|17|) 
for momentum values < k < kp, whereas for kp < k < 3kp ( s it corresponds to an 
one-electron addition exponent considered in Ref. 

While for momentum values < k < kp one reaches the same value for the exponent 
plotted in Fig. 1 in the limits m — > 0, U/t — > and U/t — > 0, m — > 0, that value is 
different for kp < k < 3kp. In this paper we always consider the limit U/t — > 0, m — > 0, 
whereas the studies of Ref. |2H1 considered the limit m — > 0, U/t — > 0. This justifies 
the different values of that exponent for kp < k < 3kp and U — given in Fig. 1 and 
in Fig. 2 of Ref. j22j, respectively, which otherwise correspond to the same exponent 
values. The U/t — > and [//£ — > oo limiting values of the exponent (|T7j) and other 
exponents obtained below are further discussed at the end of this section. 
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Figure 1. Momentum dependence of the exponents associated with the one-electron 
removal spin s = si branch line of Fig. 1 of Ref. |2H] for < k < kp and one-electron 
addition spin branch line of the same figure for kp < k < 3kp. For the one-electron 
removal case considered here that exponent is given in Eq. (|17|l . (In the figure both 
these exponents are called ( s .) We note that for U > and in the one-electron removal 
small momentum domain in the vicinity of the branch- line end point k = kp, where 
v s i(q) ~ Vsi, the exponent plotted here does not apply, since the spectral function is 
instead of the form given in Eq. (|12|l . 



Similar results are obtained for the other two one-electron removal si branch lines, 
whose exponent was not studied in Ref. We call them si,// branch lines where 

l' = ±1. Here and in the expressions provided below we use the indices t' = ±1 and 
t" = ±1 to denote contributions from processes which involve cO and si pseudofermions, 
respectively, created or annihilated at the m — > Fermi points l' 2kp = ±2/cp and 
i!' kp = ±fcp, respectively. For the s1,l' branch lines, the index l' refers to a cO 
pseudofermion particle-hole process such that a cO pseudofermion is annihilated at 
q = —L'2kp and created at q = L'2kp, where i' = ±1. The specific form of the 
general expressions for the points (k, w s i tL r(k)) belonging to the si, l' branch line in 
the m — > limit, corresponds to k = l' Akp and is given by, 

q = k + i'4:k F ; uj slill/ (q) = -e sl (q) . (18) 

The si, —1 pseudofermion singular branch line is labeled by s in Fig. 1 of Ref. [29J, 
where for uj/t < it connects the points {k = 3kp, u> = 0) and {k = 5/cf, uj = 0). In this 
case the value of the functional (jZj) is a function of k = q + l' 4kp given by the m — > 
limit of the parameter, 

2A^(q) = + I)}*; «^ = c0, si. (19) 

Use of the general expression (jl(Jj) in the m — >• limit, leads to the following expression 
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for the one-electron removal spectral function, 

#,w)«C 8li 4g)M- — - . (20) 



^ AlTy/VcO v s i ■ 

This expression corresponds to energy values just below the branch line and to bare- 
momentum values in the range —k F <q<k F and momentum values in the domains 
— 5k F < k < —3k F and 3k F < k < 5k F for l' = 1 and i! = —1, respectively. In the 
m —>■ limit, we find the following exponent expression valid for all values of U/t and 
electronic density n, 

(q)=-l+ E E {^f ~ + L ^u,si(tq° Fa u, q)Y 

av=cO, si t=±l 

1 2 

= - 1 + E {ttf- - u 'Zo + ^d),si(^k F: q)\ 

L=±l ^0 

+ E {-^7= + ^mi(^f, g)} 2 . (21) 

The dependence of the exponent (|2*T|) on the momentum is obtained by combining 
Eqs. ()18j) and (|2ip. For most of the parameter space and bare-momentum values, this 
exponent is larger than one and thus the spectral-function expression (J2*U|) does not 
describe a branch line. Consistently, for finite values of U/t the pre-factor C s \^'{q) of 
Eq. ()2U|) has smaller values than those of the pre-factor C s i{q) appearing in expression 
fllfjj) . However, for large values of U/t and bare- momentum values in the vicinity of 
l' k F such an exponent corresponds to a branch line, as it reaches values smaller than 
one. We recall that for small domains in the vicinity of the end points k = ±3k F and 
k = ±5k F the spectral function is not of the form (|2T)|) . but instead is of the general 
form given in Eq. ()12j) . 

Equivalent results are obtained for the four cO, l', l" branch lines, where l', l" = ±1. 
In this case, the specific form of the general expressions Q for the points (k, u c0: L ' t0 "(k)) 
belonging to the cO, i', l" branch lines in the m —> limit, corresponds to ko = 
i' 2k F — l" k F and reads, 

q = k + l' 2k F - l" k F ; ^do,i',i"(q) = -e c o(?) ■ (22) 

The cO, +1, +1 branch line, cO, —1, —1 branch line, and cO, —1, +1 branch line are 
labeled c, c', and c" in Fig. 1 of Ref. j2Hj, respectively, where they are represented for 
k > and uo < 0. In this case the value of the functional ((7|) is given by the m — ► 
limit of the following parameter, 

2KM = {l'^-l^-l h ^- swoNL,, q)} 2 , (23) 

where av = cO, si. 

The dependence of this quantity on the momentum k is obtained by combining 
Eqs. (I22J) and ()23|) . From use of the general expression (fTUj) in the m — > limit, we find 
the following expression for the one-electron removal spectral function, 

B(k, u) w C c o itSi //(g)(^- — j== — -) ■ 24 
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Figure 2. Momentum dependence of the exponent i125[) along the one-electron removal 
c = cO, +1, +1 branch line of Fig. 1 of Ref. (2H| for < k < kp. In the figure 
that exponent is called ( c and is given by and —3/8 for U/t — > and U/t — > oo, 
respectively. We note that for a small momentum domain in the vicinity of the branch- 
line end point k = kp, where v^(q) rs v c0 , this exponent does not apply, since the 
spectral function is of the form given in Eq. (|12|l . 



This expression corresponds to energy values just below the branch lines. In this 
expression and in the exponent expressions provided below, the bare-momentum values 
are in the range —2k F < q < 2k F . Furthermore, the corresponding momentum values 
belong to the domains —3k F < k < hp and — kp < k < 3kp for i! = i" = 1 and 
l' = l" = —1, respectively, and — 5k F < k < —k F and k F < k < 5k F for l' = —l" = 1 
and il = —l" = —1, respectively. In the m — > limit, the exponent Ceo, t', i" (?) °f 
expression (123)1 reads, 

a to pi 



C^Aq) = - 1 + E E {-"' ^ + -^f + LL " 

au=cO, si i=±l 



tOLV Si 



=±i 

1 / _ ^2 



+ {^(l + — ) + t$ s i,co(iA;jr, g)} ] . (25) 



For U/t > the exponent Ceo, +i, +i (?) of Eq. (|23j) is negative for all values of momentum, 
whereas CcO,-i,-i(?) is also m general negative, except for small values of U/t and a 
small domain of bare-momentum values. These exponents are plotted in Figs. 2 and 
3, respectively, as a function of the momentum k for k > 0, several values of U/t, 
and electronic density n = 0.59. In these figures these exponents are called ( c and 
C c ', respectively. Correspondingly, when Ceo, </, t" (?) < the weight distribution 
describes a singular branch line. 
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Figure 3. Momentum dependence of the exponent l|25(l along the one-electron removal 
charge c' = cO, —1, —1 branch line of Fig. 1 of Ref. [2H] for < k < 3kp. In the figure 
that exponent is called £c' and is given by and —3/8 for U/t — > and U/t — > oo, 
respectively. For a small momentum domain in the vicinity of the branch-line end 
point k = 3kp, where v c o(q) « v c0 , this exponent does not apply, since the spectral 
function is of the form given in Eq. (|12JI . 

In turn, the exponents Ceo, +i, -i and Ceo, -l, +i of Eq. (f2"5j) are positive. For the values 
of momentum for which these exponents are smaller than one the spectral-function 
expression describes edge branch lines. For finite values of U/t the pre-factors 
C c0i ±i i=F i(g) have in general smaller values than the pre-factors C c0i ±i i ±i(g). Moreover, 
for finite values of U/t the pre-factors C c0> +i,+i(<?) and C cQ ^i : ^i(q) are decreasing and 
and increasing functions of k, respectively, whose values are smallest for the domains 
— 3kp < k < —2kp and Ikp < k < 3kp, respectively. Again, in the vicinity of the 
branch line end points k = ±kp, k = ±3kp, and k = ±5kF, where f c o(<?) ~ ±^co, the 
spectral function is not of the form ()24|). but instead is of the general form given in Eq. 

The exponents plotted in Figs. 2 and 3 have different values in the limits m — > 0, 
U/t — > and U/t — > 0, m — > 0. This justifies the different values of these exponents 
given for U/t — > in Figs. 2 and 3 and in Fig. 3 of Ref. [22] , respectively. (Due to 
a misprint one of the U/t dependent functions contributing to the expression used in 
Ref. |2H] for the exponent given in Eq. (|23|) was replaced by its \a.rge-U/t asymptotic 
expansion; Therefore, for small finite values of U/t the values provided here for that 
exponent are also different from those given in that reference. Expression f!25|) and the 
corresponding Figs. 2 and 3 of this paper correct such a misprint. Fortunately, the only 
qualitative correction is that the exponent plotted in Fig. 3 becomes positive for small 
values of U/t and a small domain of momentum values in the vicinity of 3kp.) 

We finish this section by confirming that the momentum and energy dependence 
of the spectral-weight distribution in the vicinity of the corresponding branch lines in 
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the limits U/t —* and U/t —* oo, recovers the correct behaviors. (We recall that in 
this paper we reach the limit U/t — > by considering the limit U/t — * 0, m — > 0.) All 
expressions provided below are valid for electronic densities n such that < n < 1. 

Independently of the general exponent expressions derived by the PDT of Ref. |27j . 
we also used here the method of Refs. [221 12H| to derive the exponents associated with 
the one-electron removal spectral function expressions obtained in these references for 
U/t —>■ oo. The limiting values of the exponents obtained here fully agree with those 
obtained for U/t — > oo by use of the method of Refs. [221 E3j Thus, our general U/t 
expressions are fully consistent with the spectral function expressions found in these 
references for U/t — > oo. 

Use of the general results of Ref. [2Z| for the one-electron removal spectral function 
in the vicinity of the si branch line, which for finite values of U and energy uj is given 
by expression ([Tfijl . reveals that as U/t — > the pre-factor C s \(q) corresponding to the 
general finite-energy expression ()10|) is replaced by the weight constant of the 5-peak 
spectral function given in the second expression of Eq. |T0|) . The regimen associated 
with such a weight constant arises for very low values of U/t, where it is given by 
(1/AQ6>(?) -> 1 [27]. Its independ ence on the value of the bare-momentum q results 
from the behavior of the functional Co(g)> which for this specific branch line is such that 
Co(<?) ~^ as U/t — > for the whole corresponding domain of q values. We find that in 
the limits U/t — > and U/t—* oo, the si branch-line exponent given in Eq. (fTTjl reads, 

C.i(?) = -1, 17/* -0; C sl ( g ) = -I + 2 ( i |-) 2 , tf/f-oo, (26) 

for the q and k values of the spectral-function expression ()16|) . In turn, the pre- 
factor C s i >t corresponding to the low-energy expression (fT2|) is such that C s i >t — > 
as U/t — > 0. Thus, in that limit the second expression of Eq. (fTUj) is valid for 
the whole branch-line bare-momentum domain and the regimen associated with the 
spectral-function expression (JT2J) disappears. Moreover, according to Eq. (A2) of 
Ref. [HI], in the limit U/t — > the dispersion e s i(q) becomes the electronic spectrum 
e si(q) = — 2t[cos(g) — cos(kp)]. Consistently, according to Eq. (J2E1), the exponent (fTTjl is 
such that Csi(<?) — ^ — 1 as U/t — > for all values of q in the range < | q\ < k F and thus 
of the momentum k in the domain < | k\ < kp. Then, following the second expression 
of Eq. (jlfij) . the correct non-interacting one-electron removal spectral function is reached 
in this limit. 

For the one-electron removal si, i' branch line expression (j2T)|) . the multiplicative 
coefficient is such that C s i jt /(g) —>■ as U/t — > 0. Thus, such a branch line does not 
exist for U/t — * 0, which is the correct result. For m — > and in the limits U/t — > 
and U/t — >• oo the corresponding exponent (j2*T|) reads, 

C-M'(<Z)=3, E//t-0; UAo) = \~ l ' '^ +2 (if;) 2 ' t7*-°°>(27) 

for the g and k values of the spectral-function expression (j20l) . In the limit U/t — > oo, 
this exponent is such that Csi,t'(<?) = 5/8 for q — > i' k F , Csi,t'(9) = 3/2 for g = 0, 
and Csi, (/(<?) = 21/8 for g — > —i'kp. For the one-electron removal cO, i', t" branch line 
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expression (|24|). the multiplicative coefficient is such that Ceo, (.',(/'(?) — > as U/t — > 0, 
which again is the correct result. In the limit U/t — ► 0, the c0,l',l" branch lines 
disappear, all spectral weight being transferred over to the si branch line, which becomes 
the non- interacting one-electron removal spectrum. As the limits U/t — > and U/t — > oo 
are approached, the exponent tends to the following values, 

Cc0, t ', t »(g) = 2 - ; U/t^0; 

Cco >t ', t "(g) = - ^- + ^ ; f//t ^ oo , (28) 

for the g and values of the spectral-function expression ()24j) . Thus, in the limit 
U/t — > 0, it is given by for the branch lines such that l" = 1 and 2 for the branch 
lines such that i' l" = —1. Furthermore, for U/t — > oo the exponent is given by —3/8 
for the branch lines such that i' l" = 1 and 5/8 for the branch lines such that i! i" = — 1 . 

Hence, the one-electron removal si branch line becomes the non-interacting removal 
electronic spectrum, which corresponds to —kp < k < kp. In turn, for finite values of 
U/t the spectral weight spreads over a larger two-dimensional region of the (k, uj)- 
plane. However, most of the spectral weight is located in the vicinity of separated and 
independent cO and si branch lines and of the weak border line mentioned in Sec. IV. 
Our study provides the momentum and energy dependence of the weight distribution 
in the vicinity of such av branch lines. In the m — > limit, the maximum spread of 
the one-electron spectral- weight distribution occurs for U/t —>■ oo, where the problem 
had been already studied in Refs. [221 12H|- The U/t — > oo maximum spreading of the 
one-electron removal spectral weight at electronic density n = 1/2 is illustrated in Fig. 
1 of Ref. [22] for the spectral function B{k, uj). Our U/t — > oo expressions of Eqs. (|2l)|) 
and (J2~K|) agree with the results obtained by the method of Refs. [221 ESJ, as mentioned 
above. 

4. DISCUSSION ABOUT THE RELATION TO THE PHOTOEMISSION 
DISPERSIONS OF TTF-TCNQn AND CONCLUDING REMARKS 

An interesting realization of a quasi-lD metal is the organic charge-transfer salt TTF- 
TCNQ [2HUE1- The experimental dispersions in the electron removal spectrum of this 
quasi-lD conductor as measured by ARPES are shown in Fig. 9 (b) of Ref. 0] and Fig. 
4 of Ref. [21]. The experimental data in these figures were taken with He I radiation 
(21.2 eV) at a sample temperature of 60 K on a clean surface obtained by in situ cleavage 
of a single crystal. Instrumental energy and momentum resolution amounted to 70 meV 
and 0.07 A -1 , respectively. 

We note that the low-energy spectral properties of TTF-TCNQ involve inter-chain 
hopping and electron-phonon interactions. Thus, the ID Hubbard model PDT results 
are to be applied above the energies of these processes. The singular branch lines studied 
in Sec. Ill correspond to the (k, u;)-plane region which contains all one-electron removal 
spectral-weight singular features. In spite of the recent improvements in the resolution 
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of photoemission experiments [31 EE EI] , it is difficult to measure the exponents and the 
finest details of the electronic structure experimentally, in part due to the extrinsic 
losses that occur on very anisotropic conducting solids [33]. Thus, a crucial test for the 
suitability of the model (0) to describe real quasi-lD materials is whether the ARPES 
peak dispersions correspond to the singular branch lines and other divergent spectral 
features predicted by the PDT of Ref . [21] • 

The electronic density of TCNQ is n = 0.59 < 1. For densities in the domain 
< n < 1 and one-electron removal, the main singular spectral features predicted by 
the general PDT are of branch line type. Thus, for TCNQ the main divergent spectral 
features correspond to the singular branch lines studied in Sec. III. The only other 
singular feature is quite weak and corresponds to the lowest line of Fig. 1 of Ref. [29J . 
In the vicinity of such border line the spectral-weight distribution corresponds to the 
power-law dependence (JHJ, which is controlled by a U/t independent exponent. Due to 
its weakness such a border line does not lead to any preeminent TCNQ spectral feature. 

While the theoretical weight-distribution branch-line expressions provided in Sec. 
Ill refer to all values of U/t and n, a detailed study of the spectral-function k, u, and 
U/t dependence in the vicinity of the branch lines obtained in this paper confirms the 
validity of the preliminary predictions of Refs. [3] EH]: the electron removal spectra 
calculated for t = 0.4 eV, U = 1.96 eV (U/t = 4.90), and n = 0.59 yields an almost 
perfect agreement with the three TCNQ experimental dispersions. The exception is 
the low-energy behavior, as a result of the inter-chain hopping and electron-phonon 
interactions, as mentioned above. If accounted for a renormalization of the transfer 
integral due to a possible surface relaxation [1] , these values are in good agreement with 
estimates from other experiments [HIE]- 

The experimental TCNQ finite-energy peak dispersions of Fig. 4 of Ref. [2*H] 
correspond to the spin s = si branch line (fT4^) and charge c = cO, +1, +1 and 
d = cO, —1, —1 branch lines (J2*2*|) of Fig. 1 of that reference. Those are the main 
finite-weight singular branch lines in the one-electron removal spectral function for 
U/t = 4.90 and n = 0.59. Importantly, only these main singular features, whose line 
shape is controlled by negative exponents, lead to TCNQ peak dispersions in the real 
experiment. The exponent (fT7|) corresponds to the spin ,s = .si branch line and is 
plotted in Fig. 1 for < k < kp. The exponents (j23|) that correspond to the charge 
c = cO, +1, +1 branch line and charge d = cO, — 1, — 1 branch line are plotted in Figs. 
2 and 3, respectively. As reported in Sec. Ill, for finite values of U/t the value of 
the constant C c o, — i, -i(<z) of the spectral-function expression ([2*4^1 strongly decreases for 
momentum values such that 2kp < k < 3kp. This is consistent with the absence of 
TCNQ experimental spectral features for momentum values k > 0.597T pa 0.50 A -1 in 
Fig. 4 of Ref. [29J, along the corresponding d = cO, —1, —1 branch line of Fig. 1 of 
that reference. 

Thus, our detailed branch-line PDT analysis fully agrees with the preliminary 
theoretical results of Refs. [1] [2H] for the TCNQ problem. On the other hand, 
the theoretical predictions for the TTF dispersions presented in Ref. [2H] are very 
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preliminary. For the electronic density value corresponding to the TCNQ stacks the 
main singular spectral features are of branch-line type and the only existing border 
line is quite weak. In contrast, a careful analysis of the problem by means of the 
general PDT reveals that for the electronic density suitable to the TTF stacks the main 
singular features are both of branch-line and border-line type. Once the preliminary 
studies of TTF presented in Ref. |2H] involve the singular branch line features only, 
a very small value of U/t is predicted. However, if instead one takes into account all 
singular features provided by the PDT, the best quantitative agreement with the TTF 
experimental dispersions is reached for larger values of U/t, as confirmed elsewhere. 

In this paper we have used the exact PDT of Ref. [2Z1 to study the energy and 
momentum dependence of the one-electron removal spectral weight distribution in the 
vicinity of the singular and edge branch lines of the ID Hubbard model. A careful 
and detailed analysis of the spectral function expressions in the proximity of the charge 
and spin branch lines obtained here confirms the validity of the preliminary theoretical 
predictions of Refs. jU EH], in what the description of the band TCNQ dispersions 
observed by ARPES in the quasi- ID organic compound TTF-TCNQ is concerned. The 
TCNQ conduction band displays spectroscopic signatures of spin-charge separation on 
an energy scale of the band width. This seems to indicate that the dominant non- 
perturbative many-electron microscopic processes studied in Ref. j2Zj by means of the 
PDT and the associated scattering mechanisms investigated in Ref. [2H] control the 
unusual finite-energy spectral properties of TTF-TCNQ. The quantitative agreement for 
the whole finite-energy band width between the ID Hubbard model PDT theoretically 
predicted spectral features and the TCNQ photoemission dispersions of TTF-TCNQ 
reveals that for finite-energy the local effects of the Coulomb electronic correlations 
fully control the spectral properties of that material. Thus, we expect that the long- 
range Coulomb interactions, disorder, and impurity effects play very little role in the 
finite-energy and/or finite-temperature properties of TTF-TCNQ. That disorder and 
impurities do not play a major role is confirmed by the occurrence of spin-charge 
separation for the whole energy band width. Indeed, the presence of disorder and 
impurities would prevent the separation of the one-electron spectral-weight distribution 
in terms of spin and charge singular spectral features. 

Our present finite-energy description goes beyond the usual TLL low-energy 
investigations by means of bosonization jTT] and conformal-field theory ^7j. For low 
energy the present quantum problem is a TLL. This concept only applies to the parts 
of the one-electron spectrum of Fig. 1 of Ref. [2H] where the spectral dispersions 
can be linearized. From analysis of the figure branch lines one finds that such a 
regimen corresponds to low energies. However, our results refer to all values of the 
group velocities associated with the branch lines plotted in that figure. Thus, the spin- 
charge separation found here corresponds to the whole finite-energy band width. Only 
our finite-energy theoretical spectral features describe the experimental photoemission 
TCNQ dispersions of TTF-TCNQ, once the low-energy phase of TTF-TCNQ is not 
metallic and corresponds instead to a broken-symmetry state j2]. It follows that for the 
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present TCNQ photoemission problem, the ID physics described by the ID Hubbard 
model only becomes experimentally relevant for finite energy, where the low-energy TLL 
description does not apply. 

A detailed theoretical study of the the TTF experimental dispersions by means of 
the PDT, including consideration of both singular branch lines studied here and singular 
border lines is in progress and will be presented elsewhere. Moreover, the calculation of 
the one-electron spectral-function of the ID Hubbard model for all values of k and uo by 
use of the general PDT, which consider all contributing processes, is also in progress. 
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